Erratum: Propagation of Second sound in a superfiuid Fermi gas in the unitary limit 

[Phys. Rev. A 80, 043613 (2009)] 



Emiko Arahata and Tetsuro Nikuni 



We made some errors in Sec. IV. The correct interaction parameter is gBnB,/k-aT c = 0.807. The sound velocities 
plotted in Fig. [7] are normalized by the "Fermi velocity" defined in terms of the density and mass of bosons, i.e. 
Vf = (67r 2 nB) 1 / 3 /2Af. Due to mishandling of raw data, the values given in Fig.|8]did not correspond to our calculated 
results. The corrected version of the figure is given as Fig. [0] below. The qualitative behaviors of W\ and W2 in the 
BEC limit are quite different from those in the unitary limit, i.e. W2 > Wi in most temperatures (T > 0.4T C ). This 
change does not affect any other result nor the conclusion of the paper. 
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FIG. 0: Wi/(Wi + W 2 ) and W 2 /(Wi + W 2 ) as a function of temperature. 
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We study sound propagation in a uniform superfluid gas of Fermi atoms in the unitary limit. 
The existence of normal and superfluid components leads to appearance of two sound modes in the 
collisional regime, referred to as first and second sound. The second sound is of particular interest 
as it is a clear signal of a superfluid component. Using Landau's two-fluid hydrodynamic theory, we 
calculate hydrodynamic sound velocities and these weights in the density response function. The 
latter is used to calculate the response to a sudden modification of the external potential generating 
pulse propagation. The amplitude of a pulse which is proportional to the weight in the response 
function, is calculated the basis of the approach of Nozieres and Schmitt-Rink (NSR) for the BCS- 
BEC crossover. We show that, in a superfluid Fermi gas at unitarity, the second sound pulse is 
excited with an appreciate amplitude by density perturbations. 



I. INTRODUCTION 

Landau's two-fluid hydrodynamics describes the finite temperature dynamics of all superfluids when collisions are 
sufficiently strong to produce a state of local thermodynamic equilibrium Recent experiments have begun to 
observe sound propagation in trapped superfluid Fermi gases with a Feshbach resonance [JO]. At unitarity, the 
magnitude of the s-wave scattering length that characterizes the interactions between fermions in different hyperfinc 
states diverges {\a s \ —> oo). Owing to the strong interaction close to unitarity, the dynamics of superfluid Fermi gases 
with a Feshbach resonance at finite temperatures are expected to be described by Landau's two-fluid hydrodynamic 
equations [e| . Two-fluid hydrodynamics predicts the existence of in-phase modes in which the superfluid and normal 
fluid components move together, as well as out-of-phase modes where the two components move against to each other. 
These two sound modes in the collisional limit are referred to as first and second sounds. Of greater interest is the 
out-of-phase second sound mode, since it is a clear signal of the existence of a superfluid component. Out-of-phase 
hydrodynamic modes in strongly-interacting Fermi superfluids have been discussed theoretically in the literature. The 
propagations of first and second sound in a uniform superfluid at unitarity are discussed in Refs. 0-0]. References 



8l-lllj studied out-of-phase collective modes in trapped Fermi gases, which are more relevant to experiments. However, 
out-of-phase modes have not been observed experimentally so far. 

Experimentally, the sound wave in a highly elongated trapped gas can be excited by a sudden modification of a 
trapping potential using the focused laser beam. The resulting density perturbations propagate with a speed of sound. 



This technique was first used to probe Bogoliubov sound in a Bose-condensed gas 12| . Observed sound velocity was in 
good agreement with theoretical predictions [lsI . [l^ . Analogous sound propagations have been discussed for a normal 



Bose gas [15( . Possibility of observing propagation of first and second sound pulses in a Bose-condensed gas was also 
briefly discussed in Rcf. 15 1. Sound propagation was also studied theoretically for a normal Fermi gas in Ref. flfjl ]. 



More recently, the pulse technique was used to study sound propagation in a Fermi gas near a Feshbach resonance 
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P). In this experiment, first sound mode was observed but second sound mode was not observed. In principle, one 
should be able to probe two-fluid hydrodynamic sound modes using this technique. 

In the present paper, we discuss sound pulse propagation in a strongly interacting Fermi gas in the two-fluid 
hydrodynamic regime. In Ref. Q, the first and second sound velocities in the BCS-BEC crossover for a uniform gas 
was estimated theoretically. Reference Q also argued that both sound modes can be excited and detected both as 
density and thermal waves, but the quantitative results were not presented. Reference calculated the two-fluid 
density response spectrum in a uniform supcrfluid gas of Fermi atoms in the unitary limit, and showed that second 
sound is only weakly coupled into density response At first sight, this result seems to imply that the second 
sound cannot be excited by a density perturbation. In fact, it would not show up in Bragg scattering. However, 
we will show that second sound can be still observed by a sudden modification of the external potential generating 
pulse propagation. In this paper, we use Landau's two-fluid hydrodynamic equations to study pulse propagation in a 
unitary Fermi gas. 

In Sec. II, we discuss the linear response solutions of the Landau's two-fluid hydrodynamic equations for uniform 
superfluid gases. We show that sound pulse propagation is described in terms of the density response function. 
The amplitudes of the first and second sound pulses are explicitly expressed in terms of the weights in the density 
response spectrum. In Sec. Ill, we use the Nozieres and Schmitt-Rink (NSR) theory to calculate thermodynamic 
quantities, which are needed as inputs in our solutions of the two fluid equations. These results are used to calculate 
the temperature dependence of velocity and pulse amplitude of the second sound mode in Sec. IV. We find that second 
sound has an appreciable weight in the propagation of density pulses. For composition, in Sec. V, we calculate the 
temperature dependence of velocity and amplitude of the second sound pulse in the BEC limit using Hatrree-Fock- 
Bogoliubov-Popov (HFB-Popov) approximation. 



II. LINEAR RESPONSE SOLUTION OF LANDAU'S TWO-FLUID EQUATIONS 

In this section, we present a solution of Landau's two fluid hydrodynamic equation in the presence of external 
perturbation within the liner response theory. We review normal mode solutions. The Landau two-fluid hydrodynamic 
equations in a uniform supcrfluid are given by \v\ [ls| 

d\ 

mf t = -VP, (1) 

^+V.j = 0, (2) 

§ + V • (sv„) = 0, (3) 



dt 



m— = -V/i. (4) 



The total mass current 

mj = p s v s + PnV n , (5) 

is given in terms of the superfluid and normal fluid velocities v s and v„, as well as the superfluid and normal fluid 
densities, p s and p n . The sum of the superfluid and normal fluid densities gives the total mass density, mn = p = 
Ps + Pn- The continuity equation in Eq. @ expresses mass conservation and is always valid. Equation §3§ assumes 
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that the entropy of the fluid is carried by the normal fluid and is conserved. These equations describe reversible flow 
without any dissipation arising from transport coefficients HQ. We now consider the linearized Landau equations 
for a uniform superfluid. The linearized continuity and entropy conservation equations given by Eqs. ^fy and ([3]) are 

m ~gf + V " (flsOVs + PnQVn) = 0, (6) 

dSs , 

— + V ■ (s v n ) = 0. (7) 

Taking time derivative of Eqs. ((6|) and (O, and combining them with Eqs. ([I]) and (j4]) (in linearized forms) in 
conjunction with the thermodynamic identity nop, = SqST + 5P, we arrive at a closed set of equations in terms of the 
variables dp and Ss. Inserting the normal mode plane-wave solution 8p, 5s oc e i(q-r-"*) ; one finds ui 2 = u 2 q 2 , where u 
is given by 



u 2 = ± ^ - C 2 C 2 . (8) 

The thermodynamic quantities entering arc the adiabatic sound speed squared C 2 = (jg^J ■ f nc isothermal and the 
thermal sound speed squared C 2 = C 2 = The latter also acts as a coupling or mixing term. The 

difference between the adiabatic and isothermal sound speed squared can also be expressed as C 2 — C\ = [ ) _ . 



T 



Here, s = s/p the entropy per unit mass, and C v = T (■§§;) the specific heat per unit mass. 

We now consider the two-fluid hydrodynamics in the presence of an external time-dependent potential 5U(r,t). In 
this case, the equations for j and v s become 

= -VP - nSU, (9) 

m-^ = -VGu + SU), (10) 

Within the linear response theory, the general solution for the density fluctuation 5n can be written in terms of the 
density-density response function as 



dn(r,t) = / 7^3 / ^Xnn{^)5U(^u)e^-^\ (11) 

where SU (q, oS) — J dr J dtSU(r, t)e~ l<i ' r+ ' lult is the Fourier transform of the external potential. The density response 
function for a uniform superfluid described by the Landau's two-fluid equation is given by [Lsl ] 

Xnn(q,^) = 2 2 w 2 2^7 > ( 12 ) 

where a new velocity v is defined by 



-2 Pso dT 

v =s Q — . (13) 

PnO OS 



The two-fluid density response function in (|12[) was first derived for superfluid 4 He by Ginzburg [19j and Hohenberg 
and Martin [2^. It was first applied to weakly interacting superfluid Bosc gases by Gay and Griffin [2]]]. In the case 



of the sound propagation experiment [2j-|4|, a localized potential is applied at t > 0, while it is turned off at t = 0. 
This situation can be described as 

5U(r,t)=5U(z)6(-t). (14) 
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Here we assume that the external potential is uniform in the xy direction and is localized at z ~ 0. In this case, the 
density fluctuations at (t > 0) is given by 



MM) = ±1*1 d .SU( q)j ^e^ (t > 0), 
where x^„(q, w) = Imx nn (q, uj + irf). More explicitly, it is written as 

Imx n n(q,w + irf) 



(15) 



2 2 

ir—ZxSiuj 2 - u\q 2 ) + it—Z 2 6(uj 2 - ujq 2 ), (16) 



m m 



where 
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Using Eqs. (JTBJ) and {IS} in (fT5]l. we obtain 

<5n(z, i) = W x [5U(z - u x t) + SU{z + Ul t)} + W 2 [SU(z - u 2 t) + 5U(z + u 2 t)] . (18) 

where 

Wl= _^_ Zl= n n \i~ v i w 2 = 7 ^- 2 z 2 =- n n \ u r v \ . (19) 

2mu\ 2mu\ u\ — u 2 zraiij zraiij u\ — u 2 

The expression Eq. (|18|) describes propagation of sound pluses with the speeds u\ and u 2 with the amplitudes W\ 
and 14^ ■ We note that the above general result applies to dissipationless dynamics of all superfluid in the collisional 
hydrodynamics regime. However, details are quite different for different systems. For example, in superfluid 4 He, 
we have u 2 ~ v and hence Z 2 ~ 0. In this case only first sound can be excited by the density perturbation. In 
contrast, second sound can have an appreciable weight in the density response function in superfluid Bosc gases at 
finite temperatures. The main purpose of the present paper is to show that the second sound can be excited by the 
density perturbation of the form (|14[) in a superfluid Fermi gas at unitarity. 

The density response spectrum x"„(q, cj) in a superfluid Fermi gas at unitarity was calculated in Rcf. The 
result of Ref. showed that the weight of first sound is everywhere much larger than second sound. Second sound is 
only weakly coupled into the density response function (Z 2 ~ 0.05 is the maximum at T ~ 0.9T C ). However, as shown 
in Eq. (fT9|) . the pulse amplitude Wi involves an extra factor (1/wf) which arises due to the pulse perturbation of the 
form (|14j) . Since in general u 2 < Ux, the second sound pulse amplitude W 2 is relatively amplified, and can be much 
larger than the weight in ^"(Qi lu )- I n the following sections, we explicitly calculate u\, u 2 , and v for a superfluid 
Fermi gas using the microscopic theory. 



III. THERMODYNAMIC FUNCTIONS 

The explicit calculation of the weights Wi and W 2 in (TT9")) requires thermodynamics quantities, such as 
Pso, PnOi (dP/dp)T, s, and so on. In this section, we discuss the approximations used to evaluate these quanti- 
ties. The calculation is based on the Leggett mean-field BCS model of the BCS-BEC crossover, extended to include 
the effects of pairing fluctuations associated with the dynamics of the bound states using the approach of Nozieres and 
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Schmitt-Rink (NSR) 
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24 1 . The NSR approximation has also been used to calculate the thermodynamic properties 

In the NSR theory, the supcrfluid order 
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in the BCS-BEC crossover at both T =0 and finite temperatures 
parameter A and chemical potential p, are determined from the coupled equations Q 

m ^ 9E. 9 9c / ' 



1 



N 



El 1 



^tanh^ 
&, 2 



2 
2/3 5^ 



(20) 

Vlndetjl-^ [ S (q,^ n ) + J-ll, (21) 
£p + A 2 with £ p = e p — /i, e p = "2^-- The 



two-body s-wave scattering length is denoted as a s , and i/ n is the bosonic Matsubara frequency. The second term in 
Eq. (|21[) describes contribution from bosonic collective pair fluctuations Q, where the expression for S is given by in 
Appendix [X] The key function of interest in this paper is the thermodynamic potential, defined by 



n = -|A| 



47ra s 



i^trln [-G^+^lndet 



1 



47ra< 



G-^fc) 



ifilUr, 



Jk,k' 



(22) 
(23) 



A* i^w m + £ k 

All thermodynamic quantities of interest can be calculated once Q is given. For example, we can then calculate 
pressure by using the relation P 



■p- first. We can obtain pressure terms (dP/dp)x using numerical differentiation. 



We use the relation C 2 — C T = (§|) iL ^~ to calculate C 2 because it is difficult to calculate (dP/dp) s numerically. 
Due to the difficulty in numerical differentiation of thermodynamic quantities at low temperature T/T c < 0.2, we 
calculate u 2 , and W 2 only for T/T c > 0.2. 

The supcrfluid density p s can also be obtained from the thermodynamic potential 
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25[ . The normal fluid density 



p n associated with fcrmionic and bosonic degrees of freedom is given by the sum of their contributions: 



Pn 



2 2 dfFD(Ep) 

3m ^-^ 



E 



u 



2m d 



[l + Un° n (q,iv n ) 



an§ 2 (q,ii/ n ) 



dQ z 



[i + c/n° 2 ( q ,^„)] 



dH Q n (q, w n 



2E/n° 2 (q>„) 



an? 2 (q,w m ) 



(24) 



dQ z "™ 8Q Z 

Here, /fd(-E) = l/(e /3B + 1) is the Fcrmi-Dirac distribution function, ?7(q, iu n ) = det [1 + f/S(q, iv n )], U = and 
a supercurrent flows in the z direction with the supcrfluid velocity is v s = Q z /2m. We can obtain the supcrfluid 
density p s from the relation p s = p — p n . 

Now we present numerical results for the superfluid density p s , starting from the expression for p n given in Eq. (|24|) . 
Our calculation procedure closely follows that summarized in Refs. Q. Fig. Q] shows the calculated supcrfluid density 
p s at unitarity a s — > oo. The NSR theory does have a problem near T c near unitarity and on the BEC side of the 
crossover as a result of only considering Gaussian fluctuations. However, we consider it the best available theory for 
the thermodynamic variables at finite temperatures in the BCS-BEC crossover at the present time. 



IV. SOUND PROPAGATION 



Using the thermodynamic quantities calculated in the previous section, we discuss first and second sound propa- 
gation. Since we are interested in the unitary limit, we set l/a s — 0. In Fig. [51 we plot the sound velocities in a 
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FIG. 1: Superfluid density fraction in a uniform Fermi gas at unitarity as a function of temperature. 



uniform Fermi gas at unitarity as a function of temperature. Near the critical temperature, the second sound velocity 
approaches zero. The second sound velocity has a broad maximum around T ~ 0.9T C . The NSR-type theories devel- 
oped in Refs. [7, 22, 25, 26 1 only includes the contributions from the BCS Fermi excitations plus the bosonic pairing 



fluctuations. This leads a problem to calculate the velocity of second sound near T c 
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FIG. 2: The sound velocities in a uniform Fermi gas at unitarity as a function of temperature. Vf = hkf/m. kf is Fermi wave 
number 

In Fig. |3l we plot the temperature dependence of Wi/(Wi + W 2 ) and W 2 /(Wi + W 2 ) obtained by the NSR-type 
Gaussian fluctuation theory discussed in the previous section. One immediately sees that second sound pulse has 
an appreciable amplitude. The second sound amplitude decreases at low temperatures and increases with increasing 
temperature, before decreasing again as T c is approached. The second sound amplitude has a sharp maximum around 
T ~ 0.9T C . Figure [4] shows the perturbed density profile at T = 0.6T C for several propagation times. Since u\ > u 2l 
the second sound pulse propagates slower than the first sound pulse. We clearly see that second sound is excited 
by density perturbations in the superfluid Fermi gas at unitarity Our results show that both sound modes can be 
observed by a sudden modification of the external potential using a pulse wave. As discussed in Sec. II, although the 
weight Z 2 is very small in the unitary Fermi gas, the pulse amplitude W 2 is amplified by a factor of to because of the 
low velocity of the second sound. 

We now briefly discuss the effect of changing the scattering length. In Fig. [5j we plot W 2 /(W\ + W 2 ) as a function 
of 1/fc/a with fixing the temperature as T/T c = 0.6, where T c is the superfluid transition temperature at a given 
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FIG. 4: The perturbed density profile T = 0.6T C for several propagation times, zq = 1/kf and to = l/(fe/«j) 



1/kfa. We see that second sound pulse has an appreciable weight over a finite range in the crossover region. The 
second sound amplitude has a broad maximum around ~ 0. 
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FIG. 5: the second sound amplitude W2/CW1 + W2) as a function of 1/kfa. 
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V. THE FIRST AND SECOND SOUND IN THE BEC LIMIT 

In this section, for comparison we consider the first and second sound in the BEC limit. In the BEC limit, the 
system consists of bosonic molecules with mass M = 2m with the total number N B = N/2. The s-wave scattering 
length between molecule a B is given in terms of the atomic scattering length a s as a B — 0.62a s [271 ]. In this section, for 



simplicity we calculate the thermodynamic quantities and sound velocities of a dilute Bosc gas within the framework 
of Hatrree-Fock-Bogoliubov-Popov (HFB-Popov) approximation 28[ . This means that we assume an extreme BEC 
limit. Solving the Gross-Pitaevskii equation and the Bogoliubov equations, within HFB-Popov approximation, we 
can calculate the condensate density no and the noncondensatc density n as 



y- 

4^ v 



+ gB^o 1 I 

-gT /be (E k ) + ^ 



9b n 



- 1 



n = n B 



(25) 



where f BB (E) 



cxp(/3B)-l 



is the Bose-Einstein distribution function and 



h'k 



E, 



v e °( £ °+ 2 .9B«o) 



(26) 



Equations (j2"5")) and (|2l)|) must be calculated self-consistently. The normal fluid density is given by n n = 
^ ^2 k k 2 ^gg^ and superfluid density is given by n s = uq — n n . In Fig. [BJ we plot the temperature dependence of 




T/T c 

FIG. 6: Superfluid (condensation) density fraction in the BEC and unitary limit as a function of temperature. 



the superfluid density n s . For comparison, we also plot n s in the unitary limit. 

The thermodynamic functions can be calculated from the thermodynamic potential P = — §, where 



O = -nn Q V + -gnlV + k B T ^ In [l - e -' 3E "] 



(27) 



In Fig. [3 we plot the first and sound velocities as a function of temperature obtained from self-consistent calculation of 
Eqs. ([25|) and (|26|) . The sound velocities are normalized by the "Fermi velocity" defined in terms of the density and mass 
of bosons, i.e. v f = (6ir 2 n B ) 1/3 /2M . We fixed the parameters as gBnB/k B T c = 0.807 where T n = (jfn) 2 ^ l k B 
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within the Hartrec- 



is the BEC transition temperature. The qualitatively similar results are obtained in Ref. 
Fock approximation. In Fig. we plot the temperature dependence of W\ and Wi- When compared with Fig. [3l 
the qualitative behaviors of W\ and Wi in the BEC limit arc quite different from those in the unitary limit, i.e. 
Wi > W\ in most temperatures (T > 0.4T C ). This is mainly due to the difference of the temperature dependence in 
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FIG. 7: The first and second sounds velocities as functions of the temperature with fixed density gne/fesTc = 0.15. 
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FIG. 8: Wi/(Wi + W 2 ) and Wi/(Wi + W 2 ) as a function of temperature with fixed density g B n B /k B T c = 0.807. 



the supcrfluid density, as shown in Fig. [6j From Eq. (fT9|) . we see that the ratio W1/W2 is determined by v 2 , which is 
proportional to p s o/p n o- 



VI. CONCLUSION 



In this paper, we have discussed the propagations of the first and second sound pulses in a Fermi supcrfluid at 
unitarity. The pulse propagations are discussed in terms of the density response function obtained from Landau's two- 
fluid equations. In order to obtain all the thermodynamic quantities required for calculating the sound velocities and 
their amplitudes of the first and second sound pul ses, we use the NSR-type Gaussian fluctuation theory. The results 
for the sound velocities are consistent with Ref. ((j, . We calculated the temperature dependence of the amplitudes 
of the first and second sound mode pulses, and showed that second sound pulse has an appreciable amplitude. Our 
results show that second sound can be excited by the pulse propagation experiment and should be observed as a 
separate contribution from first sound. We hope that our results will stimulate further experiment on sound pulse 
propagation in a strongly Fermi gas in the two-fluid hydrodynamic regime. 

For composition, we also calculated the temperature dependence of velocity and amplitudes of second sound pulse 
in the BEC limit. We showed that the qualitative behaviors of Wi and W2 in the BEC limit are quite different 
from those in the unitary limit. This different is mainly due to the difference of the temperature dependence in the 
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superfluid density. 

Our work is based on a NSR-type Gaussian fluctuation theory The NSR theory does have a problem 



near T c near unitarity and on the BEC side of the crossover as a result of only considering Gaussian fluctuations. A 
more sophisticated theory will be required to obtain the results valid near T c . 
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Appendix A: Definition of H(q, iv n ) 



„, . , i / nSx+nSa + ipiSa-ngi) n^-n^ 

a(a, iv n ) = - (All 

4 V i^-nga n? 1 + n« 2 -i(iiS a -ng 1 ) 

n (l_ ^P+q/2^p-q/2 - A 2 \ E p+q/2 - Ep_ q/2 

v V E p+q / 2 E p _ q /2 J {E p+q /2 — E p _ 

x [/FD(-Ep+ g /2) - /fd(£'p-< ? /2)] 
_ L ^p+q/2^p-q/2 - A 2 \ E p+q/2 + E p _ q/2 

p \ E. p+q / 2 E p _ q /2 J {E p+q/2 + E p _ q/2 ) 2 + v 2 

x [1 - /fd(-B p+9 /2) - fME P - q /2)] , (A2) 
n o = ^p(i_ t.P+q/2tp-q/2 + A 2 \ E p+q/2 - E p _ q/2 
22 p V E p+q/2 E p _ q/2 J (E p+q/ 2 - E p _ q/2 ) 2 + vl 

x [/FD(-Ep+g/2) - /FD(-Ep-g/2)] 
_ A ^p+q/2^p-q/2 + A 2 \ + gp-g/2 

p V E p+q/2 E p _ q /2 J (E p+q/2 + E p _ q/2 ) 2 + v 2 

x [l - fFv{E p+q /2) - /fd(-E p _ 9 /2)J , (A3) 

n _ /" _ £p-g/2_\ Vn 

22 p \E p+q /2 E p _ q/2 J {E p+q/2 - E p _ q/2 ) 2 + v 2 

x [fFv{E p+q / 2 ) - fFT>{E v -q/2)\ 

_ Sp ( ^P+q/2 jp-g/2 \ Vn 

p \E p+q /2 E p „ q/2 J {E p+q/2 + E p _ q/2 ) 2 + V 2 L 

X [l - fFD{E p+q / 2 ) ~ /fd(^ p - 9 /2)J 

= (A4) 
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